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INTRODUCTION 


A  direct  ■omDar'son  ’s  "Taae  *or  'ifteen  r)asic  cases  of  oDliaue  snock-wave 

'■e’"'!  ecf:  ons  oetween  ’  nter^'erometri  c  '•esuits  otitained  at  the  Jniversity  of 

'oronto  Institute  “or  Af>rosDace  Ttuoies  'UTIAS'  10  cm  <  18  -m  Hypervel oc: ty 

Shock  "uoe  ?na  oumencal  '•esults  oDtaineo  hy  jsing  a  current  conoutatl ona  1 

nethod  “or  solving  the  Sui<»r  eauations  of  conoressible  “low.  Aoditional 

parametrized  sequences  of  calculations  are  presented  to  assess  the  utility  of 

the  present  numerical  method  in  constructing  the  various  reflection-transition 

lines  (RR  -  SMR ,  SMR  -  CMR,  CMR  -  OMR;  see  Figs.  1  &  2)  for  inviscid  flows  in  the 

shock-wave  Mach-number,  wedge-angle (M  0  )  -plane.  An  additional 

s  w 

parametrized  sequence  has  been  calculated  in  order  to  study  the  validity  of 
the  boundary-layer  displacement  theory  to  account  for  the  "von  Neumann 
paradox. " 

Over  the  past  five  '^ears,  extensive  experimental  and  analytical  data  were 
obtained  for  these  problems  (3en-Dor  4  Glass  1978,  1979,  1980;  Ando  4  Glass 
1981;  Lee  4  Glass  1982;  Shirouzu  4  Glass  1982;  Deschambault  4  Glass  1983; 
Deschambault  1984;  Hu  1984;  Hu  4  Glass  1985;  Hu  4  Shirouzu  1985;  Uheeler  and 
Glass  1985;  Wheeler  1985).  We  refer  the  reader  to  these  references  for  an 
extensive  discussion  of  the  theory  of  oblique  shock-wave  reflections,  an 
introduction  to  the  history  of  the  field,  and  further  references. 

With  the  advent  of  modern  computers,  it  has  become  possible  to  attempt 
the  computation  of  such  problems  using  finite  difference  schemes.  The  state- 
of-the-art  in  this  area  was  surveyed  in  Ben-Dor  4  Glass  1978  and  Deschambault 
4  Glass  1983;  these  authors  concluded  that  advances  in  numerical  technique 
would  be  required  before  numerical  results  could  be  viewed  with  the  same 
confidence  as  experimental  data.  The  main  object  of  this  report  is  to 
demonstrate  that  the  numerical  method  used  herein  is  sufficiently  accurate  to 
be  placed  on  a  nearly  equal  footing  with  experimental  methods  in  the  analysis 
of  perfect,  inviscid,  compressible  flows.  However,  further  development  work 
is  needed  in  the  numerical  modelling  of  nonequi 1 i hrium,  viscous  flow  fields. 

This  study  deals  exclusively  with  results  for  air  and  argon;  the 
experimental  data  for  these  results  may  be  found  in  Deschambault  1984  where 
they  are  discussed  in  detail.  Many  other  related  calculations  have  been 
performed  with  our  computer  code,  and  we  briefly  describe  them  here.  Recent 
experimental  data  for  has  also  been  obtained  and  is  reported  in  Hu  1985 


and  Hu  ^  Glass  1986.  An  analogous  numeri ca 1 -experi menta 1  study  to  the  present 
'■eoort  tor  SFg  'nay  oe  ‘ound  in  Slaz  et  al  1985.  An  interesting  problem  is 
Dosed  Dv  assuming  a  oolytropic  gas,  fixing  y.  and  and  allowing 
V  =  -atio  of  specific  neats  to  oe  a  varying  oarameter.  This  oroblem  is  well- 
suited  to  a  numerical  study  and  the  "esults  using  our  numerical  method  are 
oresented  in  lolella  ^  Glaz  1984  anc  Aerger  et  al  1985.  Finally,  a  computer 
code  has  oeen  developed  for  the  proolem  of  a  spherical  explosion  reflecting 
off  an  ideal  surface;  the  numerical  method  is  virtually  identical  to  that  used 
in  obtaining  the  results  for  this  report.  Calculations  using  this  code  are 
presented  in  Colella  h  Glaz  1984  and  Colella  et  al  1985.  The  results  of  these 
calculations  show  that  the  high  M^,  OMR  flow  fields  of  planar  oblique  shock 
wave  reflection  have  a  lot  in  common  with  the  Mach  stem  region  flow  fields  of 
the  spherical  explosion  problem  just  after  the  RR-OMR  transition,  although 
there  are  significant  structural  differences  downstream  of  the  triple  point, 
presumably  due  either  to  unsteadiness  or  the  different  boundary  conditions. 

A  portion  of  the  calculations  which  are  studied  in  this  report  have 
appeared  in  Colella  A  Glaz  1982,  1984  and  Glaz  et  al  1985.  The  latter  paper 
includes  an  expanded  discussion  of  some  of  the  overall  issues  involved  in 
comparing  experimental  results  to  approximate  solutions  of  a  perfect  inviscid 
flow.  In  this  report,  we  concentrate  on  presenting  the  complete  set  of 
calculations  including  a  discussion  of  each  comparison  or  parametric  series. 
The  plan  of  the  report  is  as  follows.  In  Section  2,  the  terminology  of 
oblique  shock-wave  reflections  is  reviewed  and  some  notation  is  defined. 
Sections  3  and  4  are  devoted  to  experimental  techniques  and  the  numerical 
method,  respectively.  In  Section  5,  the  results  are  presented  and  Section  6 
is  an  extended  summary. 


SECTION  : 


DBLIOUE  SHOCK-WAVE  REFLECTIONS 

'he  *our  types  of  pseuao-stati onary  oolique  snock-wave-'‘ef1ect1on 
oatterns  are  shown  'n  Figure  1  and  consist  of  ,'a)  regular  reflection  tRRl, 
Ibl  single  '^ach  reflection  'SHRl,  (c)  complex  ^^ach  reflection  iCMR)  ana 
;dl  dounle  ‘-lacn  reflection  ’.DWl.  Figure  1  illustrates  the  definitions  af 
wedge  angle  tr"* o  1  e-ooi nt  trajectory  angles,  x,x',  ■■'an'* ous  snock  waves.  1,  F 
R‘,  .‘•i,  H',  slip  surfaces  S,S‘  and  the  flow  regions  1-5  produced  oy  the 
foregoing  reflections,  the  angle  6  between  the  incident  I  and  reflected  R 
shock  waves  is  also  shown  as  well  as  the  angle  w'  between  R  and  the  wall  or  R 
and  the  triple-point-trajectory  angle  x-  The  bow  shock  stand-off  distances  s 
and  the  length  L,  between  the  wedge  corner  and  the  reflection  point  or  Hach 
stem  are  also  indicated.  Such  qua  ities  can  be  measured  experimentally  or 
predicted  numerically  and  provide  important  information  on  the  state  of  the 
gas  whether  frozen,  non-equilibrium  or  equilibrium  (Shirouzu  A  Glass  1982;  Hu 
1985;  Hu  i  Glass  19851. 

The  equations  of  gas  dynamics  are,  in  Cartesian  coordinates, 

2 

(pu) ^  +  (pu  +p)j^  +  (puv)^  =  0 

(1) 

(pv) ^  +  (puv)^  +  (pv^+p)^  =  0 
(pE)^  +  (puE-Hjp)^  +  (pvE+vp)^  =  0 

where  p  is  the  density,  =  (u,v)  is  the  velocity  field,  E  =  V2  (u^+v^)  +  e  is 
the  total  specific  energy,  e  is  the  specific  internal  energy,  and  p  is  the 
pressure.  The  system  is  closed  by  specifying  an  equati on-of -state  (EOS), 

p  =  P(p,e).  (2) 


p  =  (y-l)p&. 


.3) 


where  y  >  1  'S  the  ratio  of  specific  heats. 

If  ^eal  gas  and  viscous  effects  can  be  ignored  [i.e.,  equations  1),  '2'' 
hold],  the  orohlem  has  no  intrinsic  length-scale,  suggesting  the  use  of  the 
self-similar  or  pseudo-stationary  coordinate  system  (5,0)  ==  [(x-XQ)/(t-tg), 
(y"yO^/ J  where  (xg.yg)  are  the  coordinates  of  the  wedge  corner  and  tg  is 
the  time  at  which  the  incident  shock  wave  reaches  the  corner.  Following  Jones 
et  al  1951,  the  system  (1)  may  be  transformed  to  pseudo-stationary  coordinates 
and  becomes,  in  conservation  form. 


(pu)^  +  (pv)^  =  -2p 


(pu  +p)^-+  (puv)  =  -3pu 
5  n 

-2 

(puv)  +  (pv  +p)  =  -3pv 

C  n 


(4) 


where 


(puH)^  +  (pvH)^  =  -p(u^+v^)  -  2pH 


u=u-C»v  =  v-  n, 

■  1  '2  "2 
H  =  -^(u  +  V  )  +  h 


(5) 


and  h  =  e  p/p  is  the  specific  enthalpy.  We  refer  to  (u,v),  H  as  the  self 
similar  velocity  field  and  self-similar  total  enthalpy,  respectively.  In 
addition  we  define 


=  (u^  +  v^)/c^  (6) 

where  c  =  sound  speed  and  we  refer  to  m  as  the  self-similar  Mach  number.  The 
system  (4)  is,  evidently,  the  steady  Euler  equations  with  the  addition  of 
source  terms.  We  note  that  the  ratio  s/L  is  constant,  for  given  initial 
conditions,  for  self-similar  solutions  of  the  non-stati onary  equations,  just 


4 


as  s  1s  constant  tor  steady  supersonic  flow.  In  this  and  other  ways  a  change 
to  Dseudo-stati onary  coordinates  is  very  useful  in  the  analysis  of  these  ‘low 
“"ields  and  will  he  used  'n  this  study. 

In  oarticular,  the  type  of  '•eflection  oattern  ^s  a  function  of  the 
■nc‘'dent  anpcK-wave  '^ach-nunner  the  wedge  angle  a,^,  and  the  gas  ^auation 

af  state.  ~he  transition  ooundaries  in  the .  ^  ,a  -plane  “"or  ooliaue  snoc<- 

s  w 

wave  reflection  are  reproduced  from  Lee  and  Glass  (1982)  in  Figure  2  for  real 
air  and  a  polytropic  equation  of  state  with  y  =  1.40.  The  analogous  figure 
for  argon  (y  =  5/3)  may  be  found  in  this  reference.  The  construction  of  the 
transition  lines  is  based  on  various  (heuristic)  transition  criteria  and  the 
numerical  calculation  of  the  jump  conditions  at  reflection  and  triple 
points.  These  criteria,  which  have  been  the  subject  of  extensive 
investigation  in  the  literature,  are  summarized  in  Lee  and  Glass  (1984)  and 
Shirouzu  and  Glass  (1983).  In  Sec.  5,  the  numerical  results  will  be  used  to 
partly  assess  the  validity  of  some  of  these  criteria  as  well  as  the  overall 
accuracy  of  the  transition  diagram.  Figure  2. 

The  fourfold  partition  of  the  plane  illustrated  in  Figure  2  is 

quite  coarse  relative  to  the  rich  phenomenology  present  in  these  flow 
fields.  Some  other  features  that  may  be  similarly  partitioned  (see  Ben-Dor 
and  Glass  1979)  are  (a)  whether  or  not  the  reflected  shock  is  detached  or 
attached  to  the  wedge  corner;  (b)  in  the  attached  case,  whether  the  flow  at 
the  corner  is  subsonic  or  supersonic;  (c)  for  RR  whether  the  flow  is  subsonic 
or  supersonic  (in  pseudo-stationary  coordinates)  at  the  reflection  point  and 
(d)  for  SMR,  CMR  and  DMR  whether  or  not  M  "toes-out"  or  "toes-in". 

A  comprehensive  study  of  these  issues  is  beyond  the  scope  of  this  report, 
but  they  will  be  discussed  as  appropriate  in  the  comparison  of  experimental 
and  numerical  results  in  Sec.  5. 


EXPERIMENTAL  TECHNIQUES 

The  experiments  for  the  present  study  were  oerformed  in  the  UTIAS  10  cm  x 
18  cm  Hypervelocity  Shocx  Tube.  A  design,  performance  and  calibration  study 
of  the  original  facility  can  he  *"0000  in  Royer  196A.  More  recent  and  detailed 
descriptions  of  the  shocx  tube  appear  in  Rristow  1971  and  Ren-Oor  and  Whitten 
1979.  Further  details  of  the  exoeriments  associated  with  the  present  worx  can 
be  found  in  Deschambault  1984. 

3.1  EXPERIMENTAL  FACILITY. 

The  basic  shock  tube  facility  consists  of  a  1.4m  long  driver  and  a  12.2m 
channel.  The  initial  pressure  in  the  channel  can  be  easily  varied  from  near 
vacuum  to  atmospheric  conditions.  At  the  end  of  the  channel  is  a  test  section 
containing  high-quality  i nterferometric  windows  through  which  the  shock  tube 
flows  may  be  observed.  A  23-cm  diameter  fi eld-of -vi ew  Mach-Zehnder 
interferometer  (Hall  1954)  in  conjunction  with  a  giant -pulse  ruby-laser  is 
used  to  record  simultaneous  dual -wavelength  (x=694.3nm  and  347. 2nm)  infinite- 
fringe  i nterferograms  of  the  two-dimensional  flow-fields.  This  allows  the 
direct  observation  of  the  flow-field  isopycnics  (lines  of  constant  density). 
The  15ns  pulse  generated  by  the  ruby  laser  effectively  freezes  all  motion, 
thereby  producing  sharp,  clear  images. 

Two  methods  were  used  to  produce  the  incident  shock-wave  Mach-numbers  for 
the  present  study.  For  shock-wave  Mach-numbers  less  than  6  a  cold-gas  driver 
was  employed.  The  diaphragm  consisted  of  several  layers  of  my lar-polyestr- 
films.  With  the  proper  choice  of  driver  gas,  CO2  or  He,  and  diaphragm 
thickness,  the  desired  shock-wave  Mach-number  could  be  obtained  in  the  test 
gas  upon  rupture. 

For  shock-wave  Mach-numbers  greater  than  6  combustion-driver  techniques 
were  used.  Specially  scribed  stainless  steel  diaphragms  were  burst  by  the 
constant-volume  combustion  of  a  stoichiometric  mixture  of  O2  and  H2  diluted 
with  70%  He.  Combustion  was  initiated  by  the  impulsive  heating  of  a  0.38-mm 
diameter  tungsten  wire  through  the  discharge  of  a  45uF  13kV  capacitor. 

The  reflection  patterns  were  generated  by  the  impingement  of  normal  shock 
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waves  with  steel  wedges,  ’’he  wedges  were  doited  firmly  to  the  rottom  wall  of 
the  facility  to  ensure  rigidity,  fhe  sides  of  the  wedges  were  flusn  with  the 
inside  walls  and  interferometric  windows  of  the  shock-tube  test-section 

3.2  DATA  DEDUCTION  ’ECHNIQUES. 

The  infinite-fringe  interferograms  enabled  the  recording  of  small 
relative  density  cnanges  of  the  various  shock  tube  flows.  The  density 
difference  Ap  between  the  two  adjacent  fringes  of  the  same  color  is  related  to 
the  wavelength  \  of  the  light  source  (694. 3nm  and/or  347. 2nm)  and  the 
Gladstone-Dale  constant  K  (2.274  x  10"\^/kg  for  air,  \  =  589. 6nm  and 
1.574  X  10“^m^/kg  for  Ar,  X  =  694. 3nm)  and  is  expressed  by  the  relation 
Ap  =  X/KL,  where  L  is  the  depth  of  the  test  section  (10.16cm). 

To  obtain  quantitative  values  for  the  isopycnics  the  following  method  was 
employed.  From  the  initial  conditions  of  the  experiment,  i.e.,  shock-wave 
Mach -number,  wedge  angle,  initial  pressure  and  temperature,  the  thermodynamic 
states  around  the  reflection  point  for  RR  and  the  triple  point  for  MR  were 
calculated  using  two-  and  three-shock  theory  (Ando  1981).  These  were  used  as 
reference  states  from  which  all  other  density  values  could  be  obtained  using 
the  above  relation. 

The  wall-density  distribution  plots  were  obtained  directly  from  the 
interferograms.  The  origin  was  defined  to  be  the  reflection  point  of  a  RR  or 
the  foot  of  the  Mach  stem  of  a  MR.  The  corner  of  the  wedge  was  defined  to  be 
a  distance  L  from  the  origin.  All  absolute  distances  were  then  scaled  by  L 
giving  a  value  of  1  to  the  distance  from  the  origin  to  the  wedge  corner. 

Where  possible  the  center  of  the  isopycnic  intersecting  the  wedge  surface  was 
used  to  locate  the  value  of  the  density  at  that  point. 

For  some  of  the  experimental  results  presented  here,  it  was  necessary  to 
use  test  gases  with  very  low  initial  densities  and  pressures  relative  to 
atmospheric  conditions.  As  a  result,  several  interferograms  show  the  effects 
of  vibrational  nonequi 1 i brium  whi ch  must  be  taken  into  account  when  analysing 
the  corresponding  interferograms.  The  relaxation  zones  are  clearly  visible 
and  appear  as  additional  fringe  shifts  in  the  post-shock  flow-field  parallel 
to  the  frozen  incident  shock  front.  Behind  the  reflected  shock  wave,  the 
characteristic  signature  of  a  relaxing  gas  is  the  nearly  tangential  incidence 
of  the  isopycnics  and  the  reflected  shock  wave. 


SECTION  4 


NUMERICAL  METHOD 

The  numerical  results  presented  in  this  paper  have  been  calculated  -vith  a 
version  of  the  Eulerian  second-order  Godunov  scheme  for  nonstati onary  gas 
dynamics  of  a  type  considered  py  Colella  and  Woodward  1984.  The  version  of 
the  scneme  used  nere  is  presented  in  Colella  and  Glaz  1982,1983,  including  the 
modifications  required  for  non-polytropi c  gases. 

t 

The  method  is  a  finite-difference  scheme  for  systems  of  hyperholic 
conservation  laws  in  one  space-like  dimension:  for  multidimensional 
applications  such  as  the  shock -on-wedge  problem,  we  employ  operator 
splitting.  Differencing  is  in  conservation  form  and  the  numerical  fluxes  are 
computed  by  solving  zone  interface  Riemann  problems  whose  time-centred  left 
and  right  states  are  computed  from  the  characteristic  form  of  the  equations. 
This  technique  leads  to  second-order  accuracy  in  smooth  flow  and  ensures  that 
the  method  is  centred  upstream.  In  practice,  the  method  is  very  stable  and 
robust.  In  the  immediate  vicinity  of  a  strong  shock,  some  dissipation  is 
required;  this  has  been  accomplished  by  smoothly  degrading  the  scheme  to  the 
first-order  Godunov  scheme  in  such  regions.  The  degree  of  degradation  is  a 
function  of  the  shock  thickness  and  strength. 

For  argon,  we  have  used  a  perfect  (frozen)  gas  equation  of  state  with 
Y  =  5/3.  If  the  shock  tube  test  gas  was  air,  the  equation  of  state  was  chosen 

to  be  either  a  perfect  (frozen)  gas  with  y  =  7/5  or  the  Hansen  1959  real  air 

equation  of  state  as  modified  by  Deschambault  1984  for  the  present 
application.  The  efficient  solution  of  the  Riemann  problem  in  the  context  of 
our  second-order  Godunov  method  for  an  arbitrary  equation  of  state  is  treated 
in  Colella  and  Glaz  1982,1983.  Also,  these  papers  demonstrate  that  the  choice 
of  equation  of  state  has  a  substantial  influence  on  the  quantitative  numerical 
results,  as  might  be  expected. 

As  noted  in  the  preceding  section,  vibrational  non-equilibrium,  which  is 
only  temperature  dependent,  can  he  significant  for  moderate  to  high  Mach 
numbers  when  the  test  gas  is  air  (at  high  Mach  numbers  dissociation  effects 

are  also  density  dependent):  for  the  argon  cases  considered  here  we  expect  the 

gas  to  remain  frozen.  The  choice  of  an  appropriate  equation  of  state  for  the 
air  calculations  depends  mainly  on  the  vibrational  relaxation  length  1^, 
behind  the  shock  waves  I,  R,  M  of  Figure  1.  If  1^  >  1  (where  1  is  a 
characteristic  flow  length  arising  in  the  problem;  for  the  present 
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exper'inents ,  '  '  ■Llmmi,  thpn  the  gas  is  frozen  and  the  perfect  gas  equation 
of  state  is  correct.  If  1  >  1.^,  then  the  gas  is  in  equilibrium  ana  the  Hansen 
equation  of  state  for  ''eal  air  is  used.  Finally,  if  1  *  1,  then  neither  the 

frozen  nor  the  equilibrium  hypothesis  is  appropriate,  and  the  flow  is  said  to 
be  in  non-equilibrium.  We  have  numerically  treated  such  cases  as  equilibrium 
^low  “ieids  by  using  the  Hansen  eauation  of  state,  although  the  only  correct 
procedure  would  be  to  solve  an  extra  partial  differential  equation 
representing  a  rate  equation  for  vibrational  relaxation  (ana  for  dissociation 
at  high  Mach  numbers).  This  decision  will  be  an  important  issue  in  our 
discussion  of  these  cases  in  Sec.  5. 

The  computational  mesh  and  our  problem  initialization  procedure  is 
illustrated  in  Figure  3.  Note  that  these  figures  are  drawn  from  right  to  left 
to  conform  with  the  experimental  interferograms.  We  have  used  a  square 
(i.e.,  AX  =  Ay  =  constant)  mesh  for  all  ''f  the  computations  in  Sec.  5. 

Because  the  flow  is  pseudo-stationary,  the  choice  of  ax  is  immaterial. 

The  initial  data  are  taken  as  Uq,  where  U  =  (p.p.u.v,)"^  is  the  state 
vector  and  is  the  initial  shock-wave  Mach-number.  From  these  data  and  the 
given  equation  of  state,  the  post-shock  state  may  be  calculated.  To 
initialize  the  two-dimensional  calculation,  these  data  are  placed  on  the  grid 
far  upstream  (ca.  60-75  zones)  of  the  corner,  as  illustrated  in  Figure  3a; 
interpolation  of  conserved  quantities  [i.e.  L)^  =  (p  ,pu,p  v,pE  )^]  is  used  for 
zones  that  straddle  the  incident  shock.  However,  this  is  a  very  poor 
representation  of  the  numerical  shock  because  any  shock-capturing  scheme  will 
diffuse  a  shock  wave  over  two  or  more  zones  in  the  computational  mesh.  The 
resulting  structure  is  referred  to  as  a  discrete  travelling  wave  (i.e.,  a  mesh 
function  that  depends  only  on  x  -  Vt,  where  V  is  the  vector  velocity  of  the 
wave  and  equals  the  shock  speed  in  magnitude  for  a  discrete  shock  wave). 
Starting  with  any  initial  data  (e.g.,  the  one  zone  Ug  -  Uj  jump  described 
above)  satisfying  the  Rankine-Hugoniot  conditions,  the  solution  will  tend  as 
the  number  of  time-steps  becomes  large  towards  the  appropriate  discrete 
travelling  wave,  plus  other  low-amplitude  waves  that  we  refer  to  as  "starting 
error",  with  the  starting  error  separating  from  the  travelling  wave.  For  the 
present  application,  it  is  very  important  to  ensure  that  the  starting  error  is 
eliminated  before  the  shock  wave  is  allowed  to  reflect,  and  we  proceed  as 
follows.  First,  the  computer  code  is  allowed  to  run  until  the  shock  wave 


reaches  *,he  corner,  and  the  situation  in  Figure  3b  is  reached.  In  this 
‘"igure,  the  region  'mmediateiy  behind  the  shock  and  about  2-3  zones  thick  is 
the  discrete  "raveilina  wave  and  the  small  (less  than  5%)  relative  amolitucie 
disturbances  'urther  downstream  is  the  one-time  starting  error,  'he  comouter 
code  then  arbitrarily  cnanges  the  ’low  field  to  that  illustrated  in  -igure  3c, 
’.e.,  the  di'tcrete  travelling  wave  (arbitrarily  set  to  exactly  -1  zones  in  the 
computer  code';  ’s  '-etained  but  the  starting  error  is  replaced  by  the  post¬ 
shock  state  U^. 

At  this  point,  the  flow  field  becomes  truly  two-dimensional  and  the 
computer  code  is  now  run  without  further  interruption  until  the  end  of  the 
calculation  is  reached. 

The  boundary  conditions  for  this  problem,  which  are  standard,  are 
discussed  in  detail  in  Colella  and  Glaz  1983.  We  remark  here  that  our 
treatment  of  the  intersection  of  the  incident  shock  with  the  upper  or  left- 
hand  boundary  or  both  is  not  entirely  consistent  with  the  discrete  travelling 
wave  and  leads  to  the  introduction  of  a  low  relative  amplitude  (ca.  1%)  wave 
behind  the  incident  shock  at  its  intersection  with  the  boundary.  This  wave, 
which  we  call  a  boundary  error,  may  lead  to  a  rather  unaesthetic  structure  in 
the  contour  plots  and  it  can  impinge  on  the  disturbed  flow  field  behind  the 
reflected  shock.  Examples  will  be  noted  in  Sec.  5. 

All  calculations  were  performed  on  a  CRAY  I  at  Los  Alamos  National 
Laboratory,  Los  Alamos,  New  Mexico.  The  computer  code  was  designed  to  take 
significant  advantage  of  the  machine's  vector  architecture.  Each  calculation 
in  Sec.  5  required  15-40  min.  c.p.  time  with  most  in  the  range  of  20-30  min. 
Much  of  this  time  is  wasted  on  the  extra  grid  points  introduced  to  eliminate 
the  starting  error  as  well  as  grid  points  outside  the  reflected  shock.  Also  a 
fine  mesh  is  only  really  needed  in  the  Mach-stem  region.  Thus,  an  intelligent 
adaptive  mesh  structure  could  reduce  these  times  substantially. 


:ect:dn  5 

COMPUTATIONAL  RESULTS 

A  ii  rect  conparison  of  axppnmental  results  and  numerical  calculations  is 
oresentea  in  Sec.  5.1  fifteen  cases  of  oblique  shock -wave  reflections. 

-or  eignt  of  the  cases  in  air,  the  computation  has  oeen  oerforned  twice,  once 
with  a  oerr'ect  gas  EOS  with  y  =  l.A  and  once  with  the  Hansen  EOS.  "hus, 
twenty-tnree  comoutations  are  reoorted  on  in  this  oart.  In  Sec.  5.2,  the 
results  of  several  parametrized  sequences  of  calculations  are  presented  to 
demonstrate  the  capability  of  our  numerical  method  to  compute  the  correct 
transition  in  the  (M  ,9  )  -plane.  An  additional  sequence  is  presented  in  this 

S  W 

part  to  demonstrate  (upon  comparison  with  experimental  data)  the  effect  of 
boundary-layer  displacement  on  the  RR-OMR  transition. 

5.1  COMPARISON  OF  EXPERIMENT  WITH  CALCULATION. 

The  initial  conditions  for  the  fifteen  cases  are  listed  in  Table  I  along 

with  the  computational  mesh  (NX, NY)  and  the  equation  of  state  selected  for 

each  case  (and  it  is  noted  where  two  choices  of  EOS  were  made  for  a  case). 

All  four  wave  configurations  are  represented  in  the  range  of 

considered.  The  following  data  are  presented  for  each  case:  experimental 

isopycnics;  computed  isopycnics  using  the  same  density  levels  as  were  obtained 

in  the  experiment;  wall  distribution  plots,  q  vs  x/L,  with  q  =  P/P^t  p/Pq» 

e,  u  and  with  the  p/p  plot  including  a  comparison  with  experiment;  whole  flow 

0 

field  contour  plots,  using  thirty  equally  spaced  contours,  of  the  quantities 
p,e,p,M,u,v,u,v,H;  in  a  "blowup"  frame  in  the  vicinity  of  the  Mach  stem  or 
reflection  point,  contour  plots,  using  thirty  equally  spaced  contours 
of  p,e,p,M,H  are  shown  along  with  the  experimental  isopycnics,  self-similar 
streamlines,  and  a  self-similar  velocity  vector  field  plot.  For  those  cases 
involving  a  comparison  of  two  calculations  with  differing  EOS,  the  contour 
plots  of  actual  isopycnics  are  shown  together  with  the  i nterferogram,  and  an 
additional  wall  distribution  plot  is  added  comparing  p/p  vs.  x/L  for  the  two 
calculations  and  the  experiment  on  the  same  graph. 

In  order  to  assist  the  reader  in  interpreting  the  graphical  output,  we 
make  several  general  comments  here,  which  are  not  repeated  below  for  each 
case.  It  is  regretted  that  many  interesting  phenomena  are  not  commented  on  in 
the  text,  but  we  felt  it  useful  to  present  the  entire  set  of  figures.  First, 
for  those  cases  involving  two  calculations,  the  subscript  "P"  refers  to  the 


•  J»  J  t  lV  -“V  'k.*.  _'w  JnW,  ■. 


11 


oer^ect  jas  calculation  ano  the  suoscnot  "H"  refers  to  the  Hansen 
taicuiation;  -n  tne  event  tnat  sucn  a  *igure  'S  ■'eferencea  'n  the  text  ^<1  trout 
a  suDScr’ot,  the  context  determines  ^hich  'or  Oothl  ‘'■'gures  are  he^ng 
aiscusseo.  Concerning  the  contour  oiots,  the  coorainate  system  is  or'entea 
aith  the  ir'Gin  at  tne  corner  oomt,  the  <  ar  %  airectton  along  the  .veage 
surface  after  -e’' 1  ecti on ,  tne  j  or  o  airect'on  oeroenoi cu  iar  to  tne  verge  ana 
"ac'ng  jowaros  'n  the  '’gure;  However,  ve  nave  ’•everseo  ar’entation  "or  tne 
va '  1  dots  and  ■'ave  set  </L  =  1.0  at  tne  corner  ana  </L  =  1.0  at  cne 
'•erl ecti on  ooint  or  the  intersect’on  of  the  Macn  stem  vitn  the  weage 
surface.  We  regret  that  we  have  not  matched  the  length  L  in  the  oiots  of  the 
calculations  with  the  corresponding  interf erograms .  The  contour  plots  of 
those  quantities  which  may  take  on  positive  or  negative  values  use  solid 
(dashed)  lines  to  represent  positive  (negative)  contour  levels.  The  zero 
level  is  always  the  last  solid  contour.  In  particular,  the  sonic  line  in 
the  plots  can  always  be  easily  found.  The  most  important  feature  of  an 
equally  spaced  contour  plot  for  compressible  flow  is  that  discontinuities  are 
clearly  visible  because  several  contour  levels  overwrite  each  other  on  the 
plot,  at  the  location  of  the  discontinuity.  When  plotting  density  contours 
using  the  levels  prescribed  by  the  available  experimental  fringes,  this  effect 
is  still  present  but  degraded  to  varying  degrees  for  the  different  cases.  In 
particular,  density  levels  between  and  ^2*  present  at  all, 

although  in  many  cases  we  arbitrarily  inserted  extra  contours  for  aesthetic 
reasons. 

Referring  to  Figure  7  (i.e..  Case  4),  the  generic  features  of  the  various 
plots  are  discussed  (the  notation  "Figure  N"  is  used  when  it  would  not  be 
useful  to  consider  Figure  7  as  an  example).  First,  most  pages  have  a  heading 
with  certain  information:  '>S",  "ALP"  are  in  the  notation  of  the  text; 

NR  and  NZ  are  the  number  of  mesh  points  in  the  calculation  and  correspond  to 
Table  I;  "Pg"  is  p^,  the  initial  shock -tube  pressure;  KBEG  is  the  first  point 
(viewing  from  right  to  left)  in  the  x-direction  after  the  reflection  point; 
and  the  word  "PERFECT"  or  "HANSEN"  appears  to  denote  the  EOS.  Notice  that 
(NR-kbEG)  by  NZ  is  the  appropriate  aspect  ratio,  rather  than  NR  by  NZ. 

Figures  7a  and  b  are  presented  on  the  same  page  and  all  plots  are  uniformly 
labelled  according  to  the  table  appearing  along  with  these  figures.  When 
comparing  an  i nterferogram  (Figure  7a  in  this  example)  with  the  associated 
calculated  isopycnics  (Figure  7b  here),  the  effect  of  an  EOS  mismatch  in  the 
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;a'cuiaf:on  at  3^,3^  3e  s'rikinq  ana  ■^isleadinq.  Recall,  see 

-  j  ^  J 

^ect'on  t^iat;  :'’ese  se^sity  /a'ues  are  saiculated  *or  :rie  ■  nter'eroqram 
;s'nd  a  soec’^'c  s^mce  sf  eoui '  i  rr' um  lOS.  A  nume'"'cal  lonoutat'on  jssnu  a 
li'^e^ent  sro^ce  at  vDS  '  autanata  ca  1  ly  :et  aifferent  /aiues  and  smft  a:’ 
af  -.re  aoncours  away  'ron  tne’r  aor^ect  'ocacons.  An  e'<ce'''ent  examo'e  of 
:n-5  ef'ec",  'S  -'qure  '.Zbui  wne''e  nost  at  *.ne  'SODycmcs  nave  aeen  smfted  'nto 
'.’’e  "'iin’e'”  ca  ■  annex  'aye-'s  assoc'acen  ^itn  :re  ■“’■‘'ec'en  anoex  and  the  seenna 
'^acn  steri.  'ance'‘mnq  -'dures  'o  and  e,  •xe  note  :ne  ’oiiowinq:  1’ 
the  plots  in  ootn  fiqures  snow  that  the  iiistorbed  flow  is  sutasonic 
everywhere  except  in  Reqion  2  where  the  flow  is  entirely  supersonic,  (2) 
contact  or  slip  surfaces  tend  to  show  more  clearly  in  the  plots  of  e 
and  H  than  in  p  and  this  is  especially  true  of  the  boundary  of  the  vortex 
rollup  as  may  be  noted  by  comparinq  these  three  plots  in  Fiqure  7e;  of  course, 
this  effect  is  caused  by  the  different  effects  the  Ranki ne-Huqoni ot  jump 
conditions  have  on  the  number  of  contours  appearinq  in  the  shock  layers  for 
the  different  quantities  and,  it  should  also  be  noted,  the  EOS  or  the  value  of 
Y  in  the  case  of  the  polytrooic  EOS  has  a  larqe  effect,  (3)  shock  waves  can  be 
distinquished  from  slip  surfaces  by  comparinq  the  pressure  plots  with  plots  of 
3,  e,  etc.,  (^)  distinquishinq  compressions  from  rarefactions  can  lisuaHy 
be  done  with  the  pressure  contours  alone  (e.g.,  if  a  compression  steepens  to 
form  a  shock)  or  in  conjunction  with  the  wall  pressure  plot,  Fiqure  7c; 
determining  the  direction  in  which  a  wave  faces  is  often  of  interest  and  is 
not  usually  obvious  although  sometimes  the  u  plot  in  Figure  Nc  can  be  used  for 
waves  normal  to  the  wall,  (5)  the  H  plots  are  not  constant  states  because  of 
the  source  terms  in  eqn,  it  follows  from  the  Rankine-Hugoniot  conditions 
that  H  does  not  change  across  a  shock  wave  and  this  can  be  seen  clearly  in 
Figures  7d  and  e  for  the  second  Mach  stem  and  the  reflected  shock  wave, 
although  there  is  some  slight  nonmonotone  variation  inside  the  numerical  shock 
layer.  It  is  also  true  that  H  does  not  jump  across  the  incident  shock  wave 
and  the  first  Mach  stem  despite  appearances  to  the  contrary  in  figures  7d  and 
e.  Close  inspection  (not  obvious  to  the  reader  in  most  instances  in  Figures 
Nd  and  e)  reveals  that  H  is  the  same  in  each  of  Regions  0-3  at  the  first 
triple  point  and  the  variation  of  H  in  the  numerical  shock  layers  is  enough  to 
cause  this  layer  to  be  filled  in  with  several  additional  contour  levels, 
f6)  the  visual  appearance  in  Figure  7c  of  streamlines  ending  in  the  interior 
of  the  calculation  is,  of  course,  just  a  plotter  error  (the  density  of 


streanl  •!  nes  's  *i  xea  'or  '■pasons  pff  i  c’ ency  ' ,  7'  the  velocity  vector  olot, 

-igure  'e  snows  ’^ow  this  vector  jumos  across  shock  waves  and  aligns  itse’’’’ 
with  slio  surfaces  ana  the  wall  lounaary  conditions. 

"he  iisciission  of  nany  af  the  cases  -efers  to  nano  neasurenents  of 
(  and  "he  accuracy  of  these  'leasurements  's  not  usually  very  qood: 

oowever,  lifferehcas  oetween  -neasurenents  e.g.,  ’'pgaraing  the  two 
calculations  of  the  sane  case  using  different  choices  of  EOS)  is  much  more 
reliable. 

Case  1:  M  =  2.05,  3^  =  60°,  RR,  Argon.  Comparison  of  the 
experinental  and  numerical  isopycnics  (Figures  ,  p)  show  them  to  be  in  good 
agreement  with  an  error  of  about  one  fringe  at  the  start  of  the  subsonic 
region.  The  wall  density  distribution  (Figure  Ac)  disagrees  by  about  the  same 
amount.  It  nay  be  observed  that  the  density  contour  levels  curve  sharply 
towards  the  reflection  point  just  above  the  wedge  surface,  an  effect  that  ’S 
not  present  in  tne  experimental  results.  The  blowup  plots  of  M,  H  exhibit 
this  effect  as  well,  even  in  the  supersonic  region.  This  numerical  error  is 
referred  to  as  "wall  heating"  and  is  commonly  observed  in  shock  capturing 
calculations  as  shown,  for  example,  in  Noh  1976.  Wall  heating  affects  only 
the  density,  temperature,  etc.,  and  not  the  pressure  (Figure  Ad).  It  may  ne 
seen  to  account  for  part  of  the  observed  error  in  this  case,  including  the 
slight  error  in  the  value  of  the  reflected  shock  wave  density  £>2 
wall.  In  addition,  the  error  in  the  stand-off  distance  of  the  bow  shock  s, 
relative  to  the  experimental  distance  from  the  reflection  point  P  to  the 
corner  L  is  about  6.2%, 

Case  2:  m  =  1,26,  9  =  450,  RR,  Air.  Fiaures  5b  and  c  show  that 

-  s  w  ’ 

the  quantitative  agreement  between  experiment  and  calculation  is  very  poor  ‘'or 

this  case,  and  the  results  are  largely  independent  of  the  choice  of  EOS. 

Furthermore,  the  angle  between  the  wedge  and  the  reflected  shock  as  well  as 

other  gross  flow  field  quantities  are  in  substantial  error,  even  though  the 

isopycnic  patterns  are  in  excellent  qualitative  agreement.  A  possible 

explanation  for  this  severe  error  may  be  found  by  considering  the  ratio 

02/Dq  as  a  function  of  fixing  9^  =  A50.  It  turns  out  that 

o_(M  =1.26)  =  2.A9p  and  p„(M  =1.24)  =  2.09p  ,  whereas  the  calculation 
2  s  0  2  S  0 

►‘as  ^2^0^  -  2.2.  In  other  words,  the  slope  d(p2/p„)/dM^  is  so  steep  m  the 


fegion  of  interest  that  sna  ’  I  errors  :n  either  the  humerical  Tiethoo  or 


exoerl menta 1  measurenent  tan  ’eao  to  ’arge  errors  in  For  further 

liscussion  this  type  of  tons i oerat o on,  see  Hu  and  Shirouzu  1985.  Noting 

that  the  distruoed  ■"  d'm  'ieid  's  wnoily  suosomc.  Figure  5n,  ana  that  this 

tase  IS  ciose  to  the  HR-3MR  transition  ooundary,  -t  is  also  oossiole  to 

speculate  that  the  experiment  'S  actually  an  3MR  oespite  this  not  oeing 

/isible  Oh  the  •  nte'^'erogram  or  -n  the  calculations.  The  ^'u  1 1  resolution  of 

this  disagreement  should  he  possible  with  the  adaptive  mesh  version  of  our 

code,  see  Berger  et  al  1985,  used  in  a  region  of  parameter  space  around 

(M  ,3  )  =  (1.26,  45°).  A  posteriori,  it  is  seen  that  studving  the  results  for 
s  w 

both  choices  of  EOS  was  not  useful. 

Case  3;  M  =  1.50,  9  =  450,  SMR,  Air.  Comparison  of  the 

i nterferogram  and  the  calculated  isopycnics.  Figures  6a  and  b,  shows  excellent 
qualitative  agreement  and  approximately  a  one  fringe  error  quantitatively. 

This  agreement  continues  for  x/L  '  1.0  since  the  experimental  corner  flow 
field  is  inviscid.  We  measure  x  "  and  0.50  for  the  calculation  and 
experiment,  respectively.  Since  the  Mach  stem  is  only  3-4  computational  zones 
high,  possible  explanations  for  this  disagreement  include  numerical  error  due 
to  lack  of  resolution  and  the  existence  of  viscous  boundary-layer  effects  in 
the  experiment.  The  disagreement  in  the  wall  density  profiles,  Figure  6c, 
for  x/L  <  0.5  may  he  due  to  several  causes:  the  possibility  of  viscous 
effects  along  the  wedge,  the  possibility  of  the  numerical  wall-heating  error 
inte^'acting  with  the  slip  surface,  and  the  difficulty  in  precisely  locating 
the  intersection  of  a  fringe  with  the  wedge  in  the  i nterferogram.  Some  of  the 
contour  plots  in  the  blowup  frame.  Figure  5e,  illustrate  the  difficulties. 

The  different  choices  of  EOS  proved  not  to  be  important  for  this  case. 

Case  4;  M  =  3.03,  9  =  470,  OMR,  Air.  The  calculated  and 

-  s  w 

experimental  wave  patterns.  Figures  7a  and  b,  are  in  excellent  qualitative 
agreement,  including  a  relatively  sharp  slip  surface  emanating  from  the  second 
triple  point.  The  i nterferogram  shows  a  different  orientation  for  fringe  c 
and  an  extra  fringe  d  under  the  reflected  shock  between  the  two  triple  points, 
which  may  be  an  indication  that  the  gas  is  relaxing  in  this  region.  The 
differences  due  to  the  choice  of  EOS  are  small,  but  noticeable.  In 
particular,  the  values  of  x  ^nd  x*  are  close  to  the  experimental  result  for 
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Hansen  EOS  ■’ut  are  too  large  oy  anout  1'^  ^or  the  >  =  l.-l  results.  It 

snould  de  noteo,  "'owever,  tnat  d  =  -17'J  ts  very  close  to  the  RR-OMR  transition 

w 

'"'ne  ana  the  oounaary  'ayer  ■ie»‘'ect  nay  lave  had  some  e^rect  on  y  ,x  '  in  the 
exoeriment.  Also,  the  /ortex  '•oilup  is  closer  to  the  leaning  ‘-^acn  stem  for 
the  Hansen  calculation  than  ‘or  the  oer'ect  oas  ca 1 cj lat 'on ,  E’gure  ’e:  the 
■  nterH(?rogram  Joes  not  snow  the  ^olluo  moving  anead  at  all,  presumaoly  due  to 
viscous  effects.  As  is  typical  ^or  0“R  results,  the  ‘low  ^ield  is  of  mixed 
type  with  region  2  being  supersonic  and  the  remainder  being  subsonic.  Figure 
7d.  Also  typical  is  the  relative  strength  of  the  contact  surface  and  vortex 
rollup  in  the  H  plots.  Figure  7e;  the  waviness  of  this  surface  in  the 
numerical  results  is  a  hint  of  the  physical  Kel vin-Helmholtz  instability 
apparent  in  the  i nterferogram.  Concerning  the  wall  density  plots.  Figure  7e, 
the  Hansen  EOS  calculation  is  a  few  percent  high  on  the  peak  value  even  after 
correcting  the  discrepancy  in  p^/pg*  ^^6  i nterferogram,  of  course,  cannot 
exhibit  the  sharp  inviscid  peak  and  valley  in  the  rollup  region  of  the 
calculation.  The  relative  displacement  of  these  structures  between  the  two 
calculations  follows  directly  from  the  differences  in  the  calculated  values 
of  X ' •  viscous  corner  region,  as  expected,  is  not  reproduced  well  in  the 

calculations. 

Case  5:  v  =  2.65,  e  =  30°,  CMR,  Air.  Comparing  the  i nterferogram 
with  the  density  contours  in  Figures  8a  and  d,  it  may  be  seen  that  excellent 
overall  agreement  was  obtained  for  the  wave  system  except  in  the  corner 
region.  The  comparisons  using  the  expe^'imental  isopycnics.  Figures  8a  and  b, 
differ  by  a  larger  degree.  The  differences  in  the  vortex  rollup  pattern  are 
clearly  the  result  of  experimental  viscous  effects.  Figures  8d  and  e  show  a 
small  supersonic  region  at  the  triple  point,  which  is  typical  of  CMR 
results .Another  interesting  feature  of  this  flow  field  is  the  presence  of 
three  points  along  the  wedge  surface  where  u  =  0  (see  the  u  contour  plot. 
Figure  8d,  and  the  u  wall  plot.  Figure  8c);  the  first  two  occur  at  the  leading 
and  trailing  edges  of  the  vortex  rollup  pattern  and  the  third  appears  much 
further  downstream.  This  pattern  is  pervasive  (except,  see  Case  71  for  those 
Mach  reflections  with  a  vortex  rollup.  We  measure  x  "  ^nr  the  perfect 

gas  calculation  and  x  =  8.3°  for  the  Hansen  EOS  calculation  and  the 
experiment.  The  calculated  values  of  d/oq  in  the  region  0.0  s  x/L  s  0.2 
are  in  good  agreement  with  the  i nterteroaram.  Figures  Ha  and  c,  once  the 


< 


:aicuiaf:on  af  is  :orr(?ctea  "or  ^  and  the  choice  of  EOS,  and  the 

different  roHuo  oatterns  are  taken  -into  account.  The  larger  disagreements  ’n 
comDamng  the  exDerimentai  ■’soDycmcs,  rigure  3a,  and  the  wall  density 
distributions,  -’gure  3c,  downstream  of  the  vortex  nignt  oe  explaineo  by  the 
viscous  effects  oroviding  different  boundary  conditions  for  the  suosonic 
•nviscid  "low  ‘^ieid  between  the  ?'o11ud  and  the  corner  region  .wnere  these 
ef^'ects  are  suostantia  1  1 . 


Case  6:  =  5.07,  =  300,  cmR,  Argon.  The  isopycnic  patterns  are 

in  excellent  agreement,  despite  the  availability  of  relatively  few  fringes. 
Figures  9a  and  d,  except  for  the  corner  region  and  the  details  of  the  vortex 
rollup  pattern.  The  quantitative  agreement.  Figures  9a,  b  and  c  and 
measurements  of  x.  are  also  very  good  except  in  the  corner. 


Case  7-.  =  10.37,  =  10°,  CMR,  Air.  The  experimental  results. 

Figure  10a,  show  strong  '•elaxation  effects  in  the  disturbed  flow  field  behind 
the  reflected  shock  (this  is  indicated  by  the  near  tangential  incidence  of  the 
fringes  to  the  shock),  and  the  incident  shock  jump  appears  almost  in 
equilibrium.  Also,  the  wedge  surface  does  not  appear  to  be  perfectly  straight 
in  the  photograph,  which  indicates  that  the  sidewall  boundary-layer- 
diffraction  effects  may  be  significant.  There  is  reasonably  good  qualitative 
agreement  (disregarding  the  real-gas  effects)  in  the  isopycnic  patterns. 

Figure  lOd,  although  the  kink  is  more  pronounced  in  the  experiment  than  in  the 
calculation.  In  evaluating  the  wall  density  plots.  Figure  10c,  it  should  he 
noted  that  the  data  points  were  evaluated  assuming  frozen-triple-point 
conditions  while  the  calculation  implicitly  used  the  equilibrium  Hansen  EOS 
for  the  same  task.  Also,  we  estimate  x  *  13.0°  for  the  experiment  and 
measure  x  =  15.00  for  the  calculation;  the  corner  attachment  angle  is  20.5° 
for  the  experiment  and  25.5°  for  the  calculation.  The  latter  difference  is 
very  large  and  is  clearly  the  result  of  the  difference  between  an  equilibrium 
shock  jump  and  a  strongly  relaxing  shock  jump  at  the  corner.  The  former 
difference  is  probably  also  a  real  gas  effect,  and  would  have  a  strong 
influence  on  the  kink  structure.  The  vortex  rollup  patterns  are  in  remarkably 
close  qualitative  agreement,  although  we  note  the  rollup  is  closer  to  the 
leading  Mach  stem,  which  has  a  somewhat  greater  toe-out,  in  the  experiment 
than  in  the  calculation.  An  unusual  feature  of  this  flow  field  is  that  u  has 
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just  one  zero  on  the  wedge  surface,  located  at  the  leading  edge  of  vortex 
rollup.  Comparing  with  the  discussion  in  Case  5,  this  suggests  that  as 
the  -  plane  is  traversed  from  the  low  hign  '•egion  to  the  high 

,  ’  ow  9  ''egion  and  ''estricting  to  cases  *or  wmch  a  vortex  rollup  pattern 
is  oresent,  the  numoer  of  zero  crossings  of  u  along  the  wedge  surface  smoothly 
oi  furcates  oetween  one  and  three.  The  '“esults  for  Cases  14  and  15 
suostantiate  this  conjecture;  the  former  lies  near  the  transition  point  and 
has  three  zero  crossings  while  the  latter  lies  just  beyond  the  transition  (and 
the  lone  zero  crossing  ahead  of  the  vortex  rollup  is  pushed  foward  into  the 
shock  layer).  Also,  the  contact  surface  is  more  unstable  in  the  calculation 
than  in  the  experiment.  These  two  effects  are  opposite  to  those  usually 
holding  in  our  results.  Thus,  it  seems  that  quite  good  quantitative  agreement 
could  be  obtained  ^or  x/L  small  in  the  wall  density  plots  if  the  rollup 
patterns  could  be  spatially  lined  up,  the  Hansen  EOS  used  in  evaluating  the 
data,  and  the  corner  jump  conditions  changed  to  provide  the  correct  downstream 
boundary  condition  for  the  subsonic  portion  of  the  flow  field.  Figure  lOd. 

The  dip  in  Figure  10c  at  x/L  “  0.25  is  due  to  the  boundary  error. 

Case  8:  H  =  1.66,  9  =  40°,  SMR,  Air.  The  isopycnic  patterns. 

Figures  11a,  b  and  d,  as  well  as  the  wall  density  plots.  Figure  11c,  are  in 
excellent  qualitative  and  quantitative  agreement.  The  only  noticeable 
difference  between  the  two  calculations  is  that  the  value  of  o-^/oq  is  in 
better  agreement  when  using  the  Hansen  EOS  and  this  aligns  the  overall  wall 
density  plots  closer  to  the  experiment.  There  is  a  larger  error  for  x/L  large 
which  is  probably  explained  by  viscous  effects  in  the  corner  region  for  the 
experiment.  The  EOS  effect  on  the  values  of  03/00  is  worth  commenting  on  in 
detail,  since  the  small  value  of  precludes  significant  real-gas  effects. 
Assuming  a  perfect  gas  and  a  Mach  stem  normal  to  the  wedge  surface  at  the 
triple  point,  one  may  compute  03/og  =  [(y*"! )M^]/[(y -I  )M^  +  2]  where  M^  = 

M^  CSC  1^,^  and  w/2  -  (9^  +  x)  '^i’ich  implies  that  03/0^  is  sensitive  to  the 

value  of  X  3t  low  shock-wave  Mach  numbers  and/or  high  values  of  (9  •he)  and 
that  d(p-/p  )/dx  <  0.  For  this  case,  v  -  3.50  in  both  calculations  but  is 

JO 

slightly  less  in  the  Hansen  calculation  which  is  enough  to  account  for  the 
wall  density  results  presented  in  Figure  11c.  Noting  that  the  calculations 
compute  p^/p  and  account  for  the  EOS,  x,  and  any  deviation  from  normality  of 
the  Mach  stem  automatically,  while  the  experimenter  must  make  the  assumptions 


aoove  and  measure  ■<  '^and,  one  sees  chat  the  differences  in  the  various 
'■esults  are  outweighed  oy  the  agreements.  Also,  boundary  layer-displacement 
may  oe  a  factor  oecause  of  the  '“elatively  low  value  of  x* 

Case  9:  "f.  =  2.37,  0^  =  40^,  DMR,  Air.  Comparing  the  density 

contour  olots.  Figure  12d,  with  the  i nterferogram,  we  see  Chat  there  is 
excellent  overall  auaiitative  agreement  for  both  calculations.  This  agreement 
is  maintained  only  for  the  perfect  gas  calculation  when  comparing  the 
experimental  isopycnics.  Figure  12b;  of  course,  the  experimental  data 
reduction  used  a  frozen  triple  point  analysis.  The  vortex  rollup  pattern,  the 
corner  flow  field,  and  the  second  triple  point  flow  field  differ  considerably, 
however.  Taking  up  the  latter  point  first,  we  note  that  this  case  is  near  the 
CMR-OMR  transition  boundary,  irrespective  of  the  choice  of  EOS.  Also,  we 
measure  x  =  5.3°  for  the  perfect  gas  calculation,  x  =  5.00  fop  the  Hansen  EOS 
calculation  and  x  =  for  the  experiment.  Thus,  it  is  not  unreasonable  for 

the  calculations  to  contain  a  much  stronger  second  Mach  stem  and  sharper 
second  triple  point  than  the  experiment,  which  is  close  to  CMR.  The  effects 
of  boundary-layer  displacement  might  also  play  a  role.  The  experimental 
contact  surface  is  very  diffusive  and  this  effect  may  prevent  the  vortex  from 
moving  forward  towards  the  Mach  stem,  in  conjunction  with  the  presumed 
boundary-layer  effects.  The  peak  stagnation  density  (see  the  u  vs.  x/L  plots 
in  Figure  12c)  behind  the  vortex  is  substantially  higher  in  the  calculations; 
we  can  conjecture  that  this  is  due  to  the  sharper  DMR  structure  and  the 
nondiffusive  contact  surface  of  the  calculation.  In  view  of  these  effects  and 
the  differing  boundary  conditions  at  the  upstream  stagnation  point  and  the 
corner,  the  wall  density  results.  Figure  12c,  are  actually  in  good  agreement. 
The  Hansen  EOS  results  must  be  corrected  for  the  data  reduction  technique  and 
there  is  otherwise  little  difference  in  the  two  calculations.  In  particular, 
the  bunching  of  the  fringes,  in  Figure  12bf^,  at  the  second  Mach  stem  is  due  to 
the  mismatch  in  Regions  2  and  3  between  the  experiment  and  the  Hansen 
calculations. 

Case  10:  M  =  3.72,  8  =  40°,  DMR,  Air.  The  analysis  for  this  case 
follows  closely  that  for  Case  9,  although  the  i nterferogram.  Figure  13a,  is 
clearly  DMR  as  are  the  calculations.  It  is  likely  that  there  is  a  relaxation 
fringe  underneath  the  reflected  shock  between  the  two  triple  points.  We 
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^’neasure  x  ~  for  the  per’‘‘ect  gas  calculation,  x  =  o.20  ‘'or  the  Hansen  EOS 

calculation,  and  x  =  o,J0  for  the  excer’inent;  the  nifferences  in  the 
neasurenents  of  x'  similar.  Referring  to  the  contour  oiots,  -igures 
13d  ana  e,  one  sees  that  the  sonic  ’ine  ‘s  coinc‘dent  v^ith  the  second  '‘'ach 
stem;  this  always  occurs  in  our  clear  hMR  results  ana  it  is  a  useful  criterion 
in  ii sti ngui sm ng  the  CMR-OMR  transition.  Other  typical  ^low  “ield  features 
are  >1'  the  transition  of  the  second  Macn  stem  to  a  continuous  comoression 
near  its  intersection  with  the  main  contact  surface.  Figures  13d  and  e,  and 
(2)  the  existence  of  two  stagnation  points  0^  and  O2,  one  behind  the  vortex 
and  the  other  just  below  the  S-M‘  intersection  and  above  the  wedge,  see 
the  (u,v)  vector  field  plots.  Figure  13e;  note  that  the  self-similar 
streamlines  are  singular  at  these  two  points.  Also,  the  pressure  attains 
local  maxima  at  these  two  points.  Figure  13e,  and  the  u  contour  plots,  Figure 
13d,  show  u  =  0  at  Oj  and  O2.  Concerning  the  wall  density  plots.  Figure  13c, 
the  agreement  is  closer  than  it  appears  because  the  data  points  in  the 
range  0.18  <  x/L  <  0.425  need  to  be  shifted  to  the  right  to  account  for  the 
different  relative  locations  of  the  second  triple  point;  such  a  shift  lines  up 
the  plots  but  the  peaks  are  still  off  as  in  Case  9. 

Case  11;  M  =  4.62,  9  =  40®,  DMR,  Air,  The  analysis  is  similar  to 

-  s  w 

Cases  9-10.  There  is  probably  a  relaxation  fringe  underneath  the  reflected 
shock  which  is  stronger  in  the  i nterferogram.  Figure  14a,  then  for  some  of  the 
other  cases.  We  measure  x  ■  perfect  gas  calculation 

and  X  =  5.0°  for  the  Hansen  EOS  calculation  and  the  experiment;  there  are 
similar  differences  for  x'«  i’i’e  quantitative  agreement  between  the  experiment 
and  the  perfect  gas  calculation  is  very  good  away  from  the  corner  in  both  the 
isopycnic  plot.  Figure  14bp,  and  the  wall  density  plot,  Figures  14c  and  Cp.  A 
shift  of  data  points  as  in  Case  10  leads  to  nearly  exact  agreement 
for  0.16  <  x/L  <  0.35  and  the  peak  density  error  at  the  stagnation  point  is 
very  small,  relative  to  Cases  9-10.  This  is  perhaps  due  to  a  reduced  relative 
influence  of  viscous  effects  in  the  region,  although  the  i nterferogram  shows 
significant  instabilities  in  the  contact  surface  and  a  vortex  rollup  pattern 
similar  to  these  two  cases.  It  is  interesting  to  note  that  in  Cases  9-11  the 
Hansen  EOS  provides  better  agreement  with  the  experiment  in  terms  of  gross 
flow  field  features  (e.g.,  x.x')  but  worse  agreement  on  quantitative  details 
such  as  wall  density  curves  (we  are  discussing  the  situation,  of  course,  after 

t 
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the  data  nave  oeen  corrected  'or  the  cnoice  of  EOS).  The  present  case 

exemo 1 ifies  this  'act  in  that  the  Hansen  EOS  calculation  shows  exact  agreement 

on  V  but  is  oaaly  off  on  oeax  oensity  along  the  wall. 

Ease  12:  H  =  2.33,  e  =  270,  3MR,  Air.  The  agreement  oetween 

-  3  W 

calculation  and  exoeriment  is  extremely  strong  in  all  respects.  Figures  15a,  o 
and  c,  ana  is  the  Pest  of  ail  the  fifteen  cases.  Ouantitati veiy ,  the 
isopycnics  are  off  hy  about  one  fringe  and  the  wall  density  plot  shows  similar 
agreement  except  in  a  small  region  near  the  corner.  The  contact  surface 
spreads  out  in  the  experiment  and  does  not  rollup  as  much  as  in  the 
calculation. 

Case  13:  =  8.70,  8^  =  270,  Air.  The  i nterferogram.  Figure 

16a,  exhibits  substantial  real  gas  effects  and  even  the  Hansen  EOS  does  not 
model  the  isopycnic  shapes  and  locations  very  well.  The  relaxation 
length,  i  ,  is  about  0.1  x  L  for  the  incident  shock  and  the  fringes  are  at 
nearly  tangential  incidence  to  the  reflected  shock.  Also,  the  relaxing  gases 
in  the  Mach  stem  region  have  obscured  the  contact  surface  and  part  of  the 
roll-up  pattern.  The  density  contour  plot.  Figure  16d,  and  the  i nterferogram 
show  very  good  agreement.  The  rollup  patterns  substantially  agree,  although 
the  contact  surface  normal  to  the  wall  at  x/L  ■  0.02  and  the  backwards 
facing  shock  wave  Wj,  normal  to  the  wall  at  x/L  ■  0.065  in  the  calculation, 
Figure  16e,  are  either  not  resolved  or  are  lost  due  to  viscous  effects  in  the 
i nterferogram.  Both  calculation  and  experiment  exhibit  a  strong  toe-out  of 
the  first  Mach  stem;  the  kink  on  this  shock  surface  may  be  near  transition  to 
a  new  triple  point  in  view  of  the  possible  existence  of  an  extra  slip  surface 
Sg  emanating  from  this  point,  see  the  u  contour  plot,  Figure  16d,  and 
the  M  plots.  Figures  16d  and  e.  The  vector  field  plot.  Figure  16e,  shows  the 
existence  of  a  pseudo-stationary  stagnation  point  O2  near  the  intersection  of 
the  two  slipstreams,  in  addition  to  the  one  at  the  center  of  the  vortex 
indeed,  there  is  a  two-dimensional  region  around  O2  where  the  flow  appears  to 
be  stagnated.  Our  measurements  show  that  x  =  5.6°  and  x'  =  5.0°  for  the 
calculation,  and  that  x  =  7.5°  and  x’  "  7.80  for  the  experiment.  The  measured 
corner  attachment  angles  are  33.5°  and  23.0°  for  the  calculation  and 
experiment,  respectively.  This  nonequilibrium  effect  (which  apparently  is 
poorly  modelled  with  the  equilibrium  Hansen  EOS)  explains  the  large 


disagreement  near  x/L  =  i.O  in  the  wall  density  plots.  Figure  16c.  After 
correcting  the  data  for  the  Hansen  EOS,  there  are  large  errors  in  the  wall 
density  plot  in  the  range  x/L  <  0.5.  Possible  exolanations  include  the  large 
error  ■'n  downstream  doundary  condition  at  x/L  =  1.0,  the  large  difference 
’n  X  ana  x'»  viscous  effects  and  differences  in  rollup  pattern,  and  general 
'•elaxation  effects,  of  course,  the  nonequilibrium  flow  field  likely 
contributes  to  the  other  three  effects.  Overall,  real-gas  effects  have  an 
extensive  impact  on  the  flow  field  dynamics  for  this  case  and  the  equilibrium 
calculation  was  unable  to  reproduce  many  of  the  details. 

Case  14;  M  =  7.19,  e  =  200,  c/DMR,  Air.  The  i nterferogram.  Figure 
s  w 

17a,  shows  clearly  that  the  experimental  flow  field  is  neither  frozen  nor  in 
equilibrium,  including  the  disturbed  flow  beneath  the  reflected  shock.  A  more 
detailed  discussion  of  equation  of  state  and  nonequilibrium  effects  in  the 
numerical  analysis  of  this  case  is  available  in  Colella  and  Glaz  1985.  The 
triple  point  angle  x  is  nearly  in  exact  agreement,  and  the  rollup  patterns  and 
Mach  stem  toe-out  agree  qualitatively.  Figure  17d.  The  attached  shock  wave  at 
the  corner  is  bifurcated  in  the  i nterferogram  and  supersonic  in  the 
calculation,  a  possible  relaxation  effect,  although  viscous  effects  may  be 
important  too.  The  experiment  and  calculation  both  show  this  case  lying  near 
the  CMR-DMR  transition  boundary.  Lee  and  Glass  1982  conjecture  that  this 
transition  occurs  when  the  sonic  line  just  reaches  the  kink;  the  M  contours. 
Figures  17d  and  e,  bear  this  out  quite  well.  After  allowing  for  the  EOS 
correction  of  the  wall  data,  a  possible  small  shift  for  the  vortex  location, 
and  the  different  corner  structures,  the  wall  density  plots  show  surprisingly 
strong  agreement;  the  dip  at  x/L  ’  0.35  is  an  excellent  example  of  the 
computational  boundary  error. 

Case  15;  M  =  8.86,  9  =  20°,  DMR,  Air.  The  calculated  density 

contours.  Figure  18d,  and  the  interferogram.  Figure  18a,  show  good  overall 
agreement,  including  many  flow  field  details.  The  vortex  rollup  patterns  are 
very  close,  although  viscous  effects  in  the  experiment  preclude  detailed 
agreement.  The  vortex  is  pushed  toward  very  close  to  the  Mach  stem  in  both 
calculation  and  experiment;  the  calculation  shows  a  wave  interaction  W  in  this 
region  which  does  not  appear  in  any  of  the  other  cases.  This  is  seen  most 
clearly  in  the  blowup  plots.  Figure  18e.  The  details  of  this  portion  of  the 
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flow  ’"ield  arp  lost  in  the  i nterferogram  and  are  underresol  ved  ’n  the 
calculation.  However,  this  flow  field  pattern  is  reproduced  in  the 
i nterf erogram  of  Experiment  974  from  leschamhault  1984  for  wnich  = 

10.18,  9^  =  20°  and  is  '^n  air.  This  i nterferogram  is  reoroduced  here  ;n 
Figure  18f.  Vibrational  relaxation  effects  are  pervasive  in  the  experiment. 
Figure  18a,  including  the  '^acn  stem  region,  'he  failure  of  the  fringes  to 
merge  into  the  second  '^ach  stem  as  the  contours  do  in  the  calculation.  Figure 
18d,  is  probably  a  real -gas  effect.  The  corner  attachment  angle  is  27°  for 
the  calculation  and  about  21°  -  23°  for  the  experiment.  The  calculation  has  a 
supersonic  corner  and  relaxation  effects  dominate  the  experimental  results  in 
the  corner  region.  We  measure  x  “  12.20  for  the  calculation  and  x  =  ^OO  for 
the  experiment,  and  x'  =  12.50  for  the  calculation  and  11.2°  for  the 
experiment.  Overall,  nonequilibrium  effects  preclude  a  realistic  quantitative 
comparison  for  this  case. 

5.2  TRANSITION  SEQUENCES. 

Four  sets  of  parametrized  sequences  of  calculations  are  presented  in  this 
sectioh.  The  purpose  of  the  first  three  sets  of  calculations  is  to  assess  the 
potential  of  detailed  computational  results  in  constructing  oblique  shock¬ 
wave-transition  boundaries  (see  Figure  2)  and  in  validating  theories 
explaining  these  transitions.  Each  set  contains  two  sequences  of 
calculations,  one  for  a  perfect  gas  with  y  =  1.4  and  one  using  the  Hansen 
EOS.  The  following  data  is  presented  for  each  case:  whole  flow  field  contour 
plots,  using  thirty  equally  spaced  contours,  of  the  quantities  p,  in  an 
appropriate  "blowup"  frame  in  the  vicinity  of  the  triple  point  or  reflection 
point,  contour  plots,  using  thirty  equally  spaced  contours,  of  the  quantities 
p,e,p,^,H,u,  along  with  the  streamlines  and  vector  field  associated  with  the 
pseudo-stationary  velocity  (u,v).  The  purpose  of  the  fourth  set  is  to 
demonstrate  the  boundary-layer  defect  theory  by  presenting  a  parametrized 
sequence  of  inviscid  calculations  for  argon  (treated  as  a  perfect  gas  with 
Y  =  5/3)  near  the  RR-DMR  transition  boundary  and  comparing  with  an 
experimental  result.  For  this  set,  only  whole  flow  field  density  contours  are 
presented. 

Set  1 :  Here,  an  attempt  is  made  to  locate  the  SMR-CMR  and  CMR-DMR 
boundaries  for  9  =  450,  Air;  1.30  <  M  <  2.60,  perfect  gas  with  y  =  1.40; 


1.50  <  M  <  2,30,  Hansen  EOS;  in  increments  of  aM  =  0.1.  The  results  are 
5  5 

oresented  in  Figures  19  and  20,  Considering  the  olots  in  the  vicinity  of 

=  1.70,  we  see  that  the  sonic  line  has  moved  into  region  2  '"or  the  cases 

with  M  >  1.70  and  that  the  extent  of  the  supersonic  region  increases  with 

5 

increasing  shocK-wave  '^acn  nunoer.  Assuming  that  the  SMR-CMR  transition 
occurs  wnen  region  2  hecomes  suoersonic  at  the  triole  ooint  'see  Lee  and  Olass 
1982),  it  follows  that  the  =  1.70  case  is  a  CMR  and  the  cases 
where  1.30  <  <  1.70  are  SMR's  because  region  2  is  entirely  subsonic  for 

these  cases.  It  may  he  noted  that  for  =  1.30,  the  Mach  stem  M  and  the 
slipstream  S  are  only  barely  visible  and  the  case  appears  like  an  RR.  The 
differences  due  to  EOS  effect  are  not  marked  at  these  *^5  values,  but  the 
Hansen  =  1.60  results  provide  a  slightly  earlier  CMR  than  the  perfect  gas 
calculation.  The  results  agree  reasonably  well  with  the  analytic  transition 
diagram.  Figure  2.  Also,  it  would  not  be  unreasonable  for  the  reader  to  view 
Figures  19  and  20  and  take  these  transitions  at  slightly  higher  M^  values, 
which  would  have  the  effect  of  making  the  comparison  with  Figure  2  somewhat 
less  close. 

In  view  of  the  small  values  of  x  in  this  region,  it  would  be  useful  to 

restudy  these  cases  with  a  refined  mesh  in  the  triple-point  region  (using  an 

adaptive  mesh  algorithm,  Berger  et  al  1985),  thereby  substantially  eliminating 

the  effects  of  numerical  error  near  the  wall  boundary  and  allowing  sufficient 

resolution  to  separate  the  results  for  the  two  choices  of  EOS.  Also,  the 

severe  slope  de  /dM  of  the  transition  curves  at  M  •  1.70  argues  for 
w  s  s  ^ 

increased  resolution. 

We  now  consider  the  p  ,M  plots  in  the  range  2.20  <  <  2.40.  One  theory 

for  the  CMR-OMR  transition  (see  Lee  and  Glass  1982)  is  that  the  flow  at  the 

first  triple  point  should  be  supersonic  with  respect  to  the  motion  of  the 
kink.  Because  the  flow  immediately  beneath  the  reflected  shock  and  between 
the  two  triple  points  is  constant,  this  criterion  is  equivalent  to  requiring 
that  the  sonic  line  (in  pseudo-stationary  coordinates)  intersect  the  kink. 
Also,  the  sonic  line  should  have  the  same  tangent  at  the  kink  as  the  second 
^’ach  stem,  because  the  flow  is  supersonic  ahead  and  subsonic  behind  this 
discontinuity.  Finally,  the  density  contours  may  be  expected  to  begin 
coalescing  as  the  shock  wave  is  about  to  form.  Using  these  criteria,  the 

calculations  show  that  the  ^  =  2.30,  perfect  gas  case  is  a  weak  OMR  and  that 


tne  =  2.“0,  oer^'ect  gas  ;asp  is  ^  clear-cut  DMR ;  for  the  Hansen 

cal ciil at" ons,  the  h.  =  2.20  can  oe  considered  a  OVR  and  the  =  2.10  case  ’s 

a  CMR.  ’’hese  ’•esults  are  also  ’n  ''easonao’i e  agreement  with  the  analytic 

results  'or  the  oerfect-aas  transition  at  3  =  -150,  see  -ioure  2.  Note  that 

w 

r’gure  2  'noicates  that  no  OMR  can  exist  in  this  range  of  ‘"or  0^  = 

However,  ‘t  nas  oeen  *ouno  excenimenta i ly  that  the  OMR-OMR  transition  line 
meets  the  5MR-CMR  cransition  line  where  't  joins  the  RR-^^R  line.  The  exact 
shape  of  this  curve  is  not  known,  although  it  would  be  expected  to  lie  much 
closer  to  the  present  numerical  values.  Insofar  as  this  observation  is  due  to 
inviscid,  equilibrium  effects,  the  numerical  results  are  further 
corroborated.  It  would  be  of  great  interest  to  pursue  the  numerical  studies 
in  the  neighborhood  of  the  coincidence  of  the  SMR-CMR  and  CMR-DMR  lines. 

It  is  also  worth  noting  that  in  this  set,  the  isopycnic  shapes  and 
distributions  resemble  those  for  RR  until  *  1.60,  where  a  loop  exists  at 
the  wedge  corner  and  the  next  fringe  away  from  this  loop  is  bowed  towards 
it.  This  effect  becomes  increasingly  prominent  as  increases  through  the 
CMR  range,  loops  begin  to  form  near  the  slipstream  as  OMR  approaches,  and 
prominently  so  as  M^  increases  through  the  DMR  range.  For  smaller  values  of 
0^^,  such  isopycnic  distributions  can  occur  for  smaller  values  of  (see 
Figure  15,  =  2.03,  9^  =  270).  The  foregoing  gives  some  insight  into  the 

changing  overall  wave  patterns  as  the  "  plane  is  traversed. 

Set  2:  The  CMR-DMR  transition  is  studied  for  =  4.0,  Air; 

290  <  9^  <  340  perfect  gas  with  y  =  1.40;  250  <  9^  <  300,  Hansen  EOS;  in 

increments  of  A0  =  lO .  The  results  are  presented  in  Figures  21  and  22.  The 
w 

analytic  CMR-DMR  transition.  Figure  2,  for  M^  =  4.D  takes  place 

at  9  ■  320  for  a  perfect  gas  and  9  ■  260  for  the  Hansen  EOS.  The  EOS  effect 

w  w 

is  predicted  correctly,  that  is,  9^  '  320  for  a  perfect  gas  and  20°  for  the 
Hansen  EOS  so  that  this  transition  line  is  shifted  up  by  about  3°.  It  is 
worth  noting  that  the  calculated  Mach  stems  are  not  perpendicular  to  the  wedge 
at  the  triple  point;  this  is  an  assumption  in  the  analytic  calculations 
leading  to  Figure  2  (see  Lee  and  Glass  1984).  Also,  it  would  not  be 
unreasonable  to  require  the  calculated  kinks  to  clearly  sharpen  up  to  a  new 
triple  point  before  assuming  a  DMR  transition. 
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3pt  3 :  "he  ann  C‘'R-DMR  transitions  are  stuoieo  '■'or  = 

S.’S,  Air;  33  s  ^  <  103  and  22°  <  a  <  260,  oerfect  jas  witn  v  =  l.AO; 

M  W 

53  <  a  4  A3  and  153  ,  a  S  193 ,  Hansen  £OS;  in  increments  of  ae  =13.  The 

^  W  W 

'•esu’ts  are  oresenteo  'n  -iqures  23  and  2^.  "he  analytic  IMR-OfR  transit’on 

F’dure  r  'or  -  i.'5  taices  oiace  at  a  "33  for  a  oer^ect  qas  and 

at  a  ■  •'O  •■■or  the  -iansen  lOS.  Accoroino  to  our  crHena  'nvoivinc 

the  '■*  some  'me,  none  af  the  reported  calculations  with  e  ^  103  are  CMR  with 

w 

t^e  possible  exception  of  the  9  =93  Hansen  EOS  result.  Thus,  the  calculated 

w 

transitions  riiffe'*  from  the  analytic  results  by  at  least  3°.  Once  again,  none 

of  the  calculated  '^ach  stems  are  perpendicular  to  the  wedge  at  the  triple 

point.  The  analytic  C'^R-Oh'R  transition.  Figure  2,  for  =  8.75  takes  place 

at  e  ■  233  for  a  pefect  gas  and  at  0  -  I60  for  the  Hansen  EOS.  The 
w  w 

calculations  show  transition  at  9^  no  greater  than  22°  for  a  perfect  gas  and 
at  9^  ■  15-163  for  the  Hansen  EOS.  This  represents  dose  agreement.  Here  as 
well,  the  Mach  stems  are  not  perpendicular  to  the  wedge  at  the  triple  point. 

It  should  he  noted  that  the  experimental  results  are  not  in  close  agreement 
with  either  of  the  two  transition  lines  {see  Figure  2)  at  such  high  values  of 
Consequently,  new  criteria  may  have  to  be  found  so  that  better  agreement 
can  be  obtained  for  the  SMR-CW-OMR  transition  lines  (see  Deschambault  and 
Glass  1983,  and  Hu  and  Glass  1985). 

Set_£:  =  7.10,  Argon  (perfect  gas  with  y  =  5/3);  Ago  <  9^  <  553, 

in  increments  of  i9  =  13;  9  =  52.750,  53.750;  53.10  <9  <  53.50  in 

w  w  w 

increments  of  A9  =  0.1°.  The  purpose  of  this  set  of  calculations  is  to 
w 

estimate  the  inviscid  RR-DW  transition  boundary  and,  by  comparison  with 

experimental  results,  to  demonstrate  and  quantify  the  well-known  disagreement 

between  theory  and  experiment  for  this  issue  (see,  for  example,  Shirouzu  and 

Glass  1982).  An  experimental  i nterferogram  for  9^  =  Ago  and  all  of  the 

computational  results  are  presented  in  Figure  25.  Noting  the  results  in  the 

range  53.03  <  9^  <  53.50  and  comparing  with  the  experiment,  a  value  of 

Ae  ■  A.O  -  A. 53  may  be  inferred  as  the  "boundary-layer  defect"  (see  Hornung 
w 

and  Taylor  1982;  Shirouzu  and  Glass  1982;  Wheeler  and  Glass  1985)  for  the 
=  7.10  RR-DMR  transition.  We  are  referring,  in  particular,  to  the 
substantial  disagreement  concerning  the  extent  of  the  Mach  stem  region 
relative  to  the  entire  flow  field.  We  have  attempted  to  calculate  the  precise 
RR-DMR  transition  point  by  plotting  the  height  of  the  Mach  stem  relative  to  L 


against  9  for  tne  conoutati ons  and  extraDolafi  ng  tfe  curve  to  ze^o  ne'ont, 

W 

^'gure  26.  'he  result  is  9^  '  53.353  whicn  disagrees  monerate'y  w’ti^  the 
‘heoretical  '■esults  of  9  '  54.-10  ip  Log  51ass,  1982.  We  '•emarx  t^at  -.his 

er^or  niay  oe  caused  oy  an  unnoticed  nias  m  our  'neasu'"'ng  tecnmoue  done  oy 
siroiy  using  a  ruler  on  the  comouter-generated  contour  slots  of  the  olow-io 
'^acn  stem  -edion  not  snownl),  lac<  of  numerical  '■esolution  -vnpn  -.he  ‘“'aco  stem 
is  only  i-2  zones  nign,  or  a  numerical  emror  in  the  oost-snoc<  ‘low  'ielo  at 
the  wall.  In  any  case,  the  error  is  small  relative  to  the  vi scous-i nvi sci d 
difference  and  it  is  also  possible  that  the  theoretical  inviscid  roredict’on 
of  9^  '  54.40  does  not  apply  when  the  entire  disturbed  flow  field  is  taxen 
into  account.  Higher  resolution  calculations  using  an  adaptive  mesh  scheme, 
Berger  et  al  1985,  will  be  carried  out  in  an  effort  to  settle  this  issue. 

This  set  of  calculations  also  illustrates  the  dramatic  collapse  of  the  complex 
OMR-pattern  into  the  simple  RR-pattern  as  9^  changes  by  a  fraction  of  a  degree 
(see  Figures  251  and  m). 


A  conouter  :ode  las  ipen  leveloopd  ^or  ':he  ’nvisc’d,  oerfpct  las  snocK- 
■)P-wedqe  oronlen  and  the  '•esu'its  have  neen  tomoared  with  best  available 
exoe'"' -^enta  1  lata.  ’he  toae  's  'lased  an  tonteTOOrary  nethodoloqy  'n 
nune''’cal  ana'ysis  af  nyper'^ol’c  conservation  'aws,  ana  nas  only  -ecently  oee" 
a  va  1 1  a  Q 1  e . 

llood  to  excellent  qualitative  agreement  has  been  obtained  in  all  cases  o'' 
direct  comparison,  and  this  applies  to  structures  beneath  the  reflected  shocx 
such  as  the  vortex  roll-up  as  well  as  coarser  criteria  such  as  the  reflection 
oattern.  Quantitatively,  the  results  are  very  good  for  flow  fields  without 
Observable  nonequilibrium  or  viscous  effects,  except  for  Case  2.  The  error  i^ 
this  case  is  probably  a  result  of  the  relatively  large  variation  of  the 
solution  with  respect  to  small  increments  in  the  problem  parameters  in  the 
vicinity  of  the  parameter  values  defining  this  case.  When  nonequi  1  i  hr-ium  or 
viscous  effects  are  oresent,  tne  quantitative  error  can  be  10-15%  and  we  may 
recall  Case  11  which  has  a  mucn  larger,  and  unexplained  error. 

Although  not  entirely  proven,  it  appears  that  the  computer  code 
represents  a  substantial  predictive  capability  for  the  shock-on-wedge  problem 
restricted  to  inviscid,  perfect  gases.  Even  for  viscous,  real  gas  flow 
fields,  the  computational  results  provide  a  significant  amount  of  information, 
including  highly  resolved  flow-field  structures. 

Significant  non-equi  1  i brium  and  viscous  effects  have  been  demonstrated  i 
the  shock  wave  diffraction  experiments,  ^'uch  of  this  could  be  inferred 
without  the  numerical  study,  but  the  latter  can  provide  a  quantitative 
estimate  of  the  various  effects.  In  particular,  vibrational  relaxation  is 
observed  in  the  high  shock  wave  Mach  number  cases,  and  this  can  have  large- 
scale  effects  on  criteria  such  as  the  corner  attachment  angle  and  type 
subsonic  or  supersonic)  and  viscous  effects  are  important  in  determining  the 
vortex  roll-up  pattern  and  the  wedge  corner  flow  field.  Although  these 
effects  occur  in  thin  layers  or  small  regions,  they  may  have  an  effect  on  the 
quantitative  results  in  the  inviscid  portion  of  the  f ' ow  field. 

The  capability  of  the  computer  code  to  discriminate  between  very  small 

increments  in  problem  parameters  (M  ,  e  ,  and  the  equati on  of  state,  although 

S  w 

t^e  latter  has  not  been  treated  here)  has  been  demonstrated. 


By  using  Da rap’Pt ri  red  seauences  of  calculations,  't  would  Oe  ooss'iblp  to 
construct  transition  boundaries  in  the  ,  S  '-olane.  'If  rourse,  tne 

5  W 

transitions  DOtameo  would  oe  oeoendent  on  the  transition  critef'a  useo  in 
t^e'’ r  construct’ on;  our  use  of  the  sonic  criterion  in  self-sim’ar  coordinates 
snows  now  tne  'nf’mte  amount  of  data  ootentially  availanle  brom  a  calculafion 
can  Oe  ’nva’uaoie  'n  eva luati ntf  one  of  the  orooosed  criteria. 

The  discussion  of  transition  set  1  in  Section  o  illustrates  how 
parametrized  numerical  calculations  can  he  used  to  elucidate  details  of  the 
flow  field  transition  not  otherwise  available.  It  is  quite  possible  that  such 
results  will  prove  useful  in  the  discovery  of  more  precise  analytic  transition 
criteria,  in  the  future.  For  the  high  transitions,  the  inviscid  numerical 
results  provide  a  guide  for  the  analysis  of  inviscid  transition  criteria  in  a 
parameter  regime  where  analytic-experimental  agreement  has  been  relatively 
poor  and  where  nonequilibrium  phenomena  are  hard  to  avoid  in  the 
experiments.  Of  course,  the  formulation  of  transition  criteria  for  viscous, 
nonequilibrium  flow  fields  is  not  assisted  by  the  present  computer  code. 

Also  we  have  been  able  to  validate  the  conjecture  that  the  RR-Dmr. 
transition  is  offset  in  experiments  by  a  boundary-layer  defect. 

In  Section  5,  several  calculations  were  noted  where  our  analysis  could  he 
greatly  improved  with  a  more  efficient  adaptive  mesh  in  the  vicinity  of  the 
Mach  stem.  Obtaining  the  necessary  resolution  with  the  present  computer  code 
would  he  overly  expensive  if  carried  out  on  a  production  basis  for  a  large 
number  of  calculations.  Using  the  methods  of  Berger  and  Colei  la  1985,  Berger 
et  al  1985,  we  expect  to  overcome  this  problem  and  revisit  some  of  the  cases 
discussed  in  this  report.  Additionally,  we  are  working  on  techniques  to 
reduce  further  or  eliminate  the  starting  error  and  boundary  error  from  our 
results. 

In  future  work,  we  intend  to  modify  our  computer  code  and  include  an 
approximation  for  vibrational  relaxation.  We  expect  that  this  work  will 
settle  some  of  the  questions  raised  in  this  paper.  The  results  presented  here 
demonstrate,  however,  that  a  valid  approximate  solution  method  for  the  Navier- 
Stokes  equations  will  be  required  if  complete  agreement  between  experiment  and 
calculation  is  demanded.  Despite  these  shortcomings,  the  comparison  of  the 
present  numerical  simulations  with  interferometric  data  from  RR ,  SMR,  CMR  and 
OMR  experiments  are  probably  the  best  available  to  date. 
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Figure  4d.  Whole-f lowf iel d  contour-plots  -  continued. 


Figure  4.  Case  1,  M  =  2.05,  9  =  60®,  Argon,  y  =  5/3,  RR  -  Continued 
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Figure  4e  -  Blowup-frame  plots 


Figure  4.  Case  1,  =  2.05.  9^  =  60®,  Argon,  y  =  5/3,  RR  -  Continued 
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Figure  5bH.  Calculated  isopycnics  (Hansen)  using  the  experimental  fringes 
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Figure  5c.  Xall  plot  for  o/p^,  y  =  1.4  and  Hansen  calculations,  with 
experimental  data 

Figure  5.  Case  2,  =  1.26,  9  =  45’,  Air,  y  =  1.4  and  Hansen  EOS,  RR 

continued. 
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Figure  5Cp.  wall  plot  for  p/p.,  o/p^  with  experimental  data  included,  e, 

o  o 

Figure  5.  Case  2,  Mj  =  1.26,  9  =  45°,  Air,  y  -  1.4  and  Hansen  EOS,  RR  - 
continued . 
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Figure  5ep.  Blowup-frame  plots;  y  =  1.4  -  continued 


Figure  5 


Case  2,  =••  1.26,  9  =  AS",  Air,  y  =  1.4  and  Hansen  EOS,  RR 

continued. 
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Figure  5d^.  Whol e-flowfield  contour-plots;  Hansen  -  continued 


Figure  5.  Case  2,  =  1.26,  9  =  45°,  Air,  y  =  1.4  and  Hansen  EOS.  RR 

conti nued. 
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Figure  6cp. 


Figure  6 


Wall  plot  for  P/Pq,  o/dq  experimental  data  included,  e, 

y  =  1.4. 

Case  3,  =  1.50,  9^  =  45°,  Air,  y  =  1.4  and  Hansen  EOS,  SMR 

cont i nued 
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Figure  6dp.  Whol e- fl owf ield  contour- pi ots  ;  y  =  1.4  -  continued 

Figure  6.  Case  3.  =  1.50,  9.^  =  45°,  Air,  y  =  1.4  and  Hansen  EOS.  SMR 

cont i nued 
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Figure  6ep.  Blowup-frame  plots;  y  = 


Figure  6.  Case  3,  Mj  =  1.50,  9^  =  45®,  Air,  y  =  1.4  and  Hansen  EOS,  SMR 
conti nued 
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Figure  7bp.  Calculated  isopycnics  {y  =  1.4)  using  the  experimental  fringes 
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Figure  7bp.  Calculated  isopycnics  (Hansen)  using  tbe  experimental  fringes 
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-ijjr-  7.  Case  4,  Ms  =  3.03,  %  =  47°,  Air,  y  =  1.4  and  Hansen  EOS,  DMR. 
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Figure  7c-.  Wall  plot  for  p/p^.  p/Pf,  experimental  data  included,  e, 

^  u;  Y  =  1.4 

Figure  7.  Case  4,  =  3.03,  9^  =  47°,  Air,  y  =  1.4  and  Hansen  EOS,  OMR  - 

continued; 
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Figure  Ic^.  Wall  plots  for  p/Po» 
u;  Hansen. 

Figure  7.  Case  4,  =  3.03,  9^ 

continued. 


Figure  7dp.  Wbole-f lowf ipld  contour-plots;  y  = 


igure  7, 


C.3S9  4, 
cont  1  nued . 
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Figure  7ep.  Blowup-frame  plots;  y  -  1.4 


Figure  7, 


Case  4,  «  3.03.  9 

continued. 


47®,  Air,  Y  «  1.4  and  Hansen  EOS,  OMR 


Figure  7ep.  31 owup- frame  plots;  v  =  1.4  -  continued. 

Figure  7.  Case  4,  =  3.03,  -  47®,  Air,  v  =  1.4  and  Hansen  EOS,  OMR 

continued. 


Figure  7dj^.  Whole-f lowf ield  contour-plots;  Hansen 
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Case  4, 
continued. 


47“ ,  Air,  V  »  1.4  and  Hansen  EOS,  OMR 
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Figure  7d^.  Whol e- f1 owf ield  contour- pi ots ;  Hansen  -  continued. 

e  7.  Case  4,  »  3.03,  9^  =  47°,  Air,  v  *  1.4  and  Hansen  EOS,  OMR 

continued. 
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Figure  7d^.  Whol e- flowf ield  contour- pi ots  ;  Hansen  -  continued. 

Figure  7.  Case  4,  «  3.03,  9^  »  47»,  Air,  y  *  1.4  and  Hansen  EOS.  OMR 

continued. 
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Figure  7e^.  1  owup-f rane  plots;  Hanspn 


'figure  7.  Case  4.  »  3.03,  *  47*,  Air,  y  »  1.4  and  Hansen  EOS,  OMR 

cont i nued . 
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Figure  7e^.  Blowup- frame  plots;  Hansen  -  continued. 

Figure  7.  Case  4,  Mj  »  3.03,  0  «  47®,  Air,  y  »  1.4  and  Hansen  EOS,  DMR  - 
continued. 
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Figure  8bp,  Calculater*  isopycnics  (y  =  1.^)  using  the  experimental  fringes 

Figure  8.  Case  5,  =  2.65,  =  30^  Air,  y  =  1.4  and  Hansen  EOS,  CMR. 
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Figure  Bdp.  Whole-flowfield  contour- pi ots ;  y  =  -  continued. 

Figure  8.  Case  5,  =  2.65,  9^^  =  30®,  Air,  y  =  1.4  and  Hansen  EOS,  CMR 

conti nued. 
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continued. 


30°,  Air.  y 


1.4  and  Hansen  EOS 


Figure  Bd,^.  Xhol  e- FI  owf  d  contour- pi  ots  ;  iansen  -  continued. 
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Figure  9c.  Wall  plots  for  p/pg,  o/pg  with  experimental  data  included,  e, 
Figure  9.  Case  6,  *  5.07,  9^  =  30®,  Argon,  y  =  5/3,  CMR  -  continued 
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figure  lOe.  31 owup- frame  plots  -  continued.  • 

Case  ^  =  10.37,  o  -  lO” .  Air.  Hansen  EOS,  CMR  -  continued 
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Figure  11a.  Interferogram. 
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isoD/cnirs  'v  =  1.4)  using  the  e^perinental  fringes, 
* '■-  'ir  =  1.56,  n  =  40'  ,  Air,  y  =  1.4  and  Hansen  EOS,  SMR. 
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Figure  llCp.  Wall  plots  for  p/Pp.  d/pq  experimental  data  included,  e,  u 
y  =  1.4. 

Figure  11.  Case  8,  Mj  =  1.66,  9y^=  40®,  Air,  y  *  1.4  and  Hansen  EOS,  SMR  - 
continued. 
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Figure  lldp.  Whole- fiowfield  contour-plots;  y  =  1.4  -  continued. 

Figure  11.  Case  8,  =  1.66,  9  =  40®,  Air,  y  =  1.4  and  Hansen  EOS,  SMR 

continued. 
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Figure  11. 
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lld^.  Whol e- flowf ield  contour- pi ots ;  Hansen  -  continued. 

Case  8,  Mj  =  1.66,  40°,  Air,  y  *  1.^  snd  Hansen  EOS,  SMR  - 

continued.  , 
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Figure  12dp.  Whole- fiowfield  contour- pi ots ;  y  =  1.4  -  continued. 

Figure  12.  Case  9,  Mj  =  2.87,  9  =  40®,  Air,  y  =  1.4  and  Hansen  EOS,  DMR  - 

continued. 
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Figure  12dp.  Whole-flowfield  contour-plots;  y  =  1.4  -  continued. 

Figure  12.  Case  9,  =  2.87,  9^  =  40'’.  Air,  y  =  1.4  and  Hansen  EOS,  DMR 

continued. 
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Figure  12. 


Figure  12e^, 

Case  9,  Mj  = 
continued. 


Blowup- frame  plots; 
2.87,  *  40°,  Air 


Hansen  -  continued. 
y  »  1.4  and  Hansen  EOS,  OMR  - 
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Figure  13c 


p.  Wall  plots  for  p/po,  o/pg  with  experimental  data  included,  e 
Y  =  1.4. 


Figure  13.  Case  10,  =  3.72,  9  =  40“,  Air,  y  =  1.4  and  Hansen  EOS,  OMR 

continued. 
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^i^ure  13'1_|.  //hoi  e- ovitfield  contour- pi  ot$  ;  Hansen  -  continued 

13.  Case  10,  =  3.72.  =  40%  Air,  v  =  1.4  and  Hansen  EOS, 

conti nued. 
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Figure  13d^.  Whole-flowfield  contour- pi ots  ;  Hansen  -  continued. 

Figure  13,  Case  10,  M  =  3.72,  9  =  40“,  Air,  y  =  1.4  and  Hansen  EOS,  DMR 

continued. 
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//all  plots  for  o/Pg,  o/og  experimental  data  included,  e 

>  =  1 . 

:.jse  11.  =  4.62,  *4.^  =  do®.  Air,  y  =  1.4  and  Hansen  EOS,  OMR 

c jnt 1 nued . 
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Figure  14clp.  Whol e- fl owf iel d  contour-plots;  y  =  continued. 

Figure  14.  Case  11.  =  4.62,  9^  =  40°,  Air,  y  =  1.4  and  Hansen  EOS,  OMR 

continued. 
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Figure  14ep.  Blowup- frame  plots;  y  = 

Figure  14.  Case  11,  =  4.62,  9^  =  40°,  Air,  y  =  1.4  and  Hansen  EOS,  DHR 

continued. 
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Figure  14du.  VJhole-flowfield  contour-plots;  Hansen. 


Figure  14.  Case  11,  Mj  =  4.62,  9^  =  40®,  Air,  y  =  1.4  and  Hansen  EOS,  DMR 
continued. 
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Figure  14e^.  B1 owup- frame  plots;  Hansen. 

Case  11,  =  4.62,  9^^  =  40”,  Air,  y  =  1.4  and  Hansen  EOS,  OMR 

continued . 
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Figure  17a. 

Interferogram. 

Figure  17b.  Calculated  isopycnics  using  the  experimental  fringes. 


Figure  17.  Case  14,  =  7.19,  0^  =  20°,  Air,  Hansen  EOS,  C/DMR. 
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Figure  17d.  Whole-flowfield  contour-plots  -  continued. 


Figure  17.  Case  14,  =  7.19,  ^  =  20®,  Air,  Hansen  EOS,  C/DMR  -  continued 
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Figure  17e.  Blowup- frame  plots. 

Figure  17.  Case  14,  »  7.19,  9^  »  20“ ,  Air,  Hansen  EOS,  C/DMR  -  continued, 
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Figure  17e.  Blowup- frame  plots  -  continued. 


Figure  17.  Case  14,  »  7.19,  0^  «  20”,  Air,  Hansen  EOS,  C/DMR  -  continued. 
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Figure  18d.  Whole-flowfield  contour-plots. 

Figure  18.  Case  15,  =  8.86,  ^  =  20",  Air,  fiansen  EOS,  DMR  -  continued 
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Figure  19.  Transition  set  1,  8^,  =  45°,  y  *  1.^  -  continued 
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Figure  19.  Transition  set  1.  ^  *  45®,  y  =  1.4  -  continued 
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Figure  19.10a.  =  2.20,  whole-fiowfield  contour-plots. 

Figure  19.  Transition  set  1,  =  45®,  y  »  1.4  -  continued 
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Figure  19.12a.  M.  =  2.40,  whole-flowfield  contour- plots. 


Figure  19.  Transition  set  1,  9,^  =  45°,  y  .  1.4  -  continued. 
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Figure  19.  Transition  set  1.  =  45®,  y  =  1.4  -  continued 
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2.60,  blowup- frame  plots  -  continued 
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Figure  19.  Transition  set  1,  ^  =  45°,  y  = 
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1.4  '  continued 
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Figure  20.1b.  =  1.50,  blowup-frame  plots. 

Figure  20.  Transition  set  1.  ^  =  45°,  Hansen  -  continued. 
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Figure  20.2b.  =  1.60,  blowup-frame  plots. 

Figure  20.  Transition  set  1,  9,^  =  45®,  Hansen  -  continued. 
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wup- frame  plots. 

5°,  Hansen  -  continued. 
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Figure  20.  Transition  set 
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Figure  21.3a.  8^  «  31®,  whole-flowfield  contour- plots. 
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Figure  21.  Transition  set  2,  M  »  4.0,  y  -  1.4  -  continued. 
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Figure  21.  Transition  set  2,  Mj 


4,0,  Y  *  -  continued 
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Figure  21.6a.  9^  =  34®,  whole-flowfield  contour-plots. 

Figure  21.  Transition  set  2,  =  4.0,  y  =  1.4  -  continued 
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Figure  22.1b.  ^  »  25®,  blowup- frame  plots  -  continued. 

Figure  22.  Transition  set  2,  M.  =  4.0,  Hansen  -  continued 
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Figure  23.3b.  9^  =  8”,  blowup- frame  plots  -  continued. 


Figure  23.  Transition  set  3,  M.  =  8.75,  y  =  1.4  -  continued. 


*.*«•»**'.  b s  y 


:  9.30 

:u=33r 

IMILAR  ENThRLR’' 

'30  _EvE 

1.98eM0  r0  6.15C*10  STEP  5.09E*0e  LflOELS  *I.O£-Oe 

STRFfiMLiNES  0F  SELF-SIMILAR  FL0W 


-3..1«*OS  T0  6.29E*(M  STEP  I.Z6E*01  lSBElS  *1.0£-:3 

self-similar  VEL0CITY  VECT0RS 


V 

V 

V 

1  1 

1  i  i 

i  i 

4  * 

; 

; 

i  i 

1  i  i 

*  * 

i  • 

\ 

v 

; 

; 

V  l 

i  i  t 

\ 

i  i 

;  i  i 

4  4 

4  4 

V 

i  i 

i  1  1 

i  4 

4  4 

\ 

V 

\  i 

i  i  i 

1  i 

4  4 

\ 

\ 

\ 

i  ( 

i  i  i 

4  4 

4  4 

\ 

\ 

\ 

4  4 

i  t  4 

.  , 

4  4 

Figure  23.6b.  *  22° ,  blowup- frame  plots. 

Figure  23.  Transition  set  3,  M.  =  8.75,  v  =  1.4  -  continued 
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Figure  24.2a.  =  6°,  whol e- fl owf ield  contour- pi ots . 


Figure  24.  Transition  set  3,  M  =  8.75,  Hansen  -  continued 
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Figure  24.9a.  9^,^  =  18® ,  whole-flowfield  contour- plots. 


Figure  24.  Transition  set  3,  =  8.75.  Hansen  -  continued. 
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Figure  25.  Transitiorr  set  4,  Mj  *  7.10,  y  »  5/3,  density  contour  plots  - 
continued. 
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Figure  26.  piot  of  DMR  Mach  stem  height  versus  e  ,  extrapolated  to  zero  height 
for  RR(h/L  *  0  for  =  53,850),  h/L  *  0  for  540  is  a  numerical  result 
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